Ising and Spin orders in Iron-based Superconductors 
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Motivated by recent neutron scattering experiments, we study the ordering of spins in the iron- 
based superconductors La(Oi_xFx)FeAs, assuming them in proximity to a Mott insulator in the 
phase diagram. The ground state of the parent system with a; = is a spin density wave with 
ordering wave vector Q = (0, tt) or (7r,0). Upon raising the temperature, we find the system to 
restore SU(2) symmetry, while an Ising symmetry remains broken, explaining the experimentally 
observed lattice distortion to a monoclinic crystal structure. Upon further temperature increase, 
the spins finally disorder at a second transition. The phase transition driven by doping with charge 
carriers similarly splits into an 0(3) transition, and an Ising transition with ^ = 3 at larger doping. 



After more than two decades of prevailing in condensed 
matter physics, copper-based high temperature super- 
conductors have very recently given in to their iron-based 
cousins [l|, y, y, y, Isl, la, [3] • The newly discovered mate- 
rials M(Oi_xFx)FeAs with M being rare earths such as 
La, Sm, have similar layered structure with stacked FeAs 
planes, sandwiched with La and O. Transport measure- 
ments show that the ground state of the undoped parent 
system is not an insulator, and LDA calculations have 
identified both a srnall electron pocket and a hole pocket 
at the Fermi level [8|. However, it has been argued that 
the system is actually close to a Mott insulator, and a 
lot of physics can be studied in a similar manner as the 
copper-based high Tc family, especially in the undoped 
system [9|. In our current work, we will study the mag- 
netism of these materials. Although the true unit cell 
of the FeAs plane contains two Fe ions, because of the 
staggered out-of-plane distribution of As ions (Fig. [1]), 
we are only interested in the magnetic Fe ions and so will 
use a unit cell with one Fe ion, unless stated otherwise. 

Recent neutron scattering experiments have shown 
that by lowering temperature the undoped material first 
undergoes a structural phase transition at 150K, with a 
distortion from tetragonal structure to monoclinic struc- 
ture, followed by a spin ordering phase transition at 134K 
developing stripe order at (tt, 0) [lO| . The observed lat- 
tice distortion and spin density wave (SDW) pattern are 
depicted in Fig. [TJ In the superconducting material with 
X = 0.08, the SDW order is not observed and, surpris- 
ingly, the distortion is absent as well, which suggests that 
the lattice distortion is driven by the development of the 
spin order. Interestingly, however, the SDW order at 
low doping only appears at a lower temperature (134K) 
than the lattice distortion (150K). We argue here that 
the SDW ordering is preceded, both in temperature and 
doping, by the breaking of an Ising symmetry in the ef- 
fective spin model, and that this is responsible for the 
observed lattice distortion. 

Ref. [Qj] argued that the undoped material is described 
by either a5' = lorS' = 2 spin model with nearest and 
next-nearest neighbor couplings Ji , J2 that depend on 
the competition between the onsite Hubbard interaction 



and the Hunds rule: 



H — 2_^ JiSi ■ Sj + 2_^ JiSi ■ Sj 



(1) 



<hj> 



«ij> 



Upon doping, this J1-J2 model has ^2,2. 

dxy superconductivity as J2 increases into the regime 



y-2+idxy and then 



where the insulator has (tt, 0) SDW order [ll|. There 
is also a much weaker interlayer coupling Jj_, which is 
necessary to stabilize the spin order. It was suggested 
by first principle calculations that both Ji and J2 are 
large and antiferromagnetic I4I, and Ref. 13| showed 
that J2 ~ 2Ji. It is well-known that when Ji < 2J2, the 
classical ground state manifold of model H]) is S'^ S'^, 
because the two sublattices of the square lattice will each 
form a Neel order (ni and 7x2), and the ground state en- 
ergy is independent of the relative angle between these 
two Neel vectors. However, quantum or thermal fluc- 
tuations lift the degeneracy, leading to parallel or an- 
tiparallel alignment of the two sublattice Neel vectors 
IJ, |l5|, ll6| . If we define 0(3) vectors (j)i as fii with soft- 
ened unit-length constraint, the long wave-length field 
theory reads: 



a— 1 fj.—x,y 

L' = "t(j)idxdy ■ (j)2 - a{(j)i ■ $2f ■ 



(2) 



In the above equation we have absorbed the overall en- 
ergy scale into (pa ■ The parameter r is tuned by tempera- 
ture, 7 ~ J\l Ji, a has contributions from both quantum 
and thermal fluctuations: a ~ J1/J2 x i^Q + ITT/J2), 
coefficients 7q and 7t are given in Ref. [15| . L' contains 
all sublattice couplings preserving the square lattice sym- 
metry. The latter rules out the term pi ■ pi, but allows 
for the couphng (0i ■ (pi)'^ ■ 

The ground state manifold of the field theory ([2]) is 
S^ (8) Z2, and the Z2 order can be described by the 
Hubbard-Stratonovich field $, which couples to (pi ■ (p2' 
L' = -^{$1 ■ $2) + $^/(4a). The ordered state with 
<i> = 1 ($ = -1) corresponds to the (tt, 0) ((0,7r)) 
SDW order. States with Ising $ order, but only short- 



range SDW order first appeared [la, llll in the quan- 
tum theory of H for 5* = 1, 2. If the coupling a is 
relevant, an Ising variable a can be introduced directly 
as 4>2 = (j4>i- Ref. [15| showed that thermal fluctu- 
ations renormalize the anisotropy mixing 7 to zero at 
long wavelength, so that at large scales the Lagrangian 
([2]) can also be viewed as the low energy field theory 
of the following Ising-0(3) model on the square lattice: 
H = X]<i 7> J{^+o'iO'j)ni 'TT-j- The 0(3) vector n denotes 
either of ni or 7x2, and the coarse-grained mode of a is 
precisely the Ising field $ introduced before. The easy- 
plane version of the Ising-0(3) model, dubbed the Ising- 
XY model, has been used widely as an effective model 
for the fully frustrated XY model on the s quar e lattice 
and the triangular lattice [H [H S [lH, [H, [H, [H . 



Note that the Ising order $ 7^ docs not imply 0(3) 
order; however, because of the system is invariant un- 
der exchanging 0i and (j)2, an 0(3) order in (^1^2 implies 
Ising order. Therefore the transition temperature of the 
Ising order is necessarily higher than that of the 0(3) 
symmetry breaking. If we consider a purely two dimen- 
sional system, at finite temperature there is only a 2d 
Ising transition separating an Ising ordered phase and a 
disordered phase since a uniform 0(3) order cannot ex- 
ist at finite temperature in dimensions smaller than 3. 
The transition temperature can be estimated roughly as 
Tci/(J20i) ~ C^/a^, £, is the correlation length of the 2d 
0(3) order at the transition, and ^ ja^ is a factor gained 
from integrating out the 0(3) order parameters. A more 
precise estimate of the Ising transition temperature for 
the J\ — J2 nrodel can be found in Ref. 15| , with T^x given 

by T = 0.13^ X ^^, in the large-S" limit. The Ising 
order breaks the 7r/2 rotation symmetry of the square lat- 
tice, indeed, an order parameter $ = 1 implies that the 
spins tend to be aligned parallel along x but antiparallel 
along y. This Ising order favors a lattice contraction in 
y direction, i.e., towards the orthorhombic structure in 
Fig. [TJ The lattice distortion thus exists even in the ab- 
sence of a uniform 0(3) order, but it necessarily requires 
the Ising order. A similar mechanism was proposed for 
the lattice distortion in the cuprates [25| . 

The interlayer coupling J^ will drive the 2d Ising tran- 
sition to a 3d Ising transition, but since it is much weaker 
than the intralayer couplings, it will not move the tran- 
sition temperature significantly. However, the interlayer 
coupling stabilizes an 0(3) ordered phase at finite tem- 
perature, assuming the interlayer coupling is small, the 
transition temperature can be estimated as follows: The 
correlation length of the 2d 0(3) nonlinear sigma model 
scales as ^/a ~ exp(27r/3), with g ^ j ^a , ^na is the mag- 
netic moment of the SDW order observed at zero temper- 
ature in units of the Bohr magneton, which is empirically 
found to be only m ~ 0.36; the interlayer coupling J± 
grows under renormalization, and becomes unperturba- 
tive when J±/{J2m'^) x C^/a^ ~ 1, which will lead to the 
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FIG. 1: (a), the lattice structure at room temperature. The 
gray circles are Fe ions, the green and red circles are As 
ions above and below the Fe plane, respectively. The dashed 
square is the two-Fe unit cell. (6), the lattice structure be- 
tween 150K and 134K after the Ising order is developed: the 
thick lines represent the bonds between antiparallel aligned 
spins, but no uniform spin order is formed. The one-Fe unit 
cell is orthorhombic, but the two-Fe unit cell has a 3d mono- 
clinic structure, as was seen in Ref. [10|. (c), the (0, tt) spin 
order below 134K. (d), the global phase diagram as a function 
of temperature and doping x. The blue curve represents the 
Ising transition, the red curve represents the 0(3) transition. 



transition temperature rc2 ^ \-n J2rr? j \a{vr? J2I J \}) [2d |. 
J2 is evaluated to be ^^lOOOK in Ref. |133, using the tran- 
sition temperature from Ref. [10], the interlayer coupling 
Jj_ is estimated to be of order lO^^J2m^, which can be 
neglected as compared with other interactions, unless we 
are very close to a critical point. The small value of the 
moment m is probably due to quantum fluctuations at 
zero temperature, since the system can be close to quan- 
tum phase transitions. Close to, but above Tc, the corre- 
lation length of the system scales like in the 3d 0(3) uni- 
versality class, but once £,z/a ~ [{T — Tc2) / J±]^'' shrinks 
to 1, the system crosses over to two dimensional critical 
behavior. The fact that the lattice distortion observed 
in experiments [10] occurs at a temperature which is rel- 
atively small compared to the exchange interaction J2 
[13| is probably due to the proximity to quantum phase 
transitions, which is consistent with the small magnetic 
moment observed at low temperature [101 ]. The phase 
diagram is depicted in Fig. [1] Notice that the lattice 
is compressed along the ferromagnetic order direction of 
the SDW, as observed experimentally J27| . 

Quantum phase transitions - In LaOi_xFxFeAs, the 
SDW order vanishes as a small amount of extra carriers 
are introduced by doping, meanwhile the superconduc- 



tor state emerges, implying the presence of one or more 
quantum phase transitions as a function of doping. A 
tentative quantum critical point in these systems has al- 
ready been studied experimentally in a series of samples 
SmOi_xFxFeAs |28|. Since the nature of the supercon- 
ductor is not yet clear, however, in the present work we 
focus on the quantum phase transitions of the spin system 
discarding the presence of superconductivity. In terms of 
the itinerant fermions, the SDW at (tt, 0) can be under- 
stood from the large susceptibility arising from the loca- 
tion of electron and hole pockets in the Brillouin zone: 
there are low energy electron-hole pair excitations at the 
(tt, 0) wavevector (Fig[2|). As extra electrons are doped 
into this system, these low energy excitations disappear 
rapidly because of the unequal sizes of the electron and 
hole pockets. The SDW order parameter cannot decay 
into a particle-hole pair excitations preserving both mo- 
mentum and energy, because the SDW wave vector (tt, 0) 
does not connect two pairs of points on the Fermi surface 
for finite x (Fig. [2|). After integrating out electrons we 
would obtain the following z = \ Lagrangian: 

2 

i—1 fi—T,x,y 
L' = -f$id,dy$2+ll\$i\^\$2\^-a{$i-$2f, (3) 

which contains no dissipative term. The first three terms 
of the Lagrangian describe the two copies of 3d 0(3) 
systems on the two sublattices. The first term in L' mixes 
(pi and (p2, and its scaling dimension is 



A[7] = D - (2 + D - 2 + r/) = -r/. 



(4) 



ry — 0.0375 [29| is the anomalous dimension of (f) at the 
3d 0(3) universality class, therefore the 7 term is irrel- 
evant. The second term in L' is allowed by symmetry 
and hence will be generated under renormalization. Its 
scaling dimension can be evaluated as 

A[7i] =D- 2A[|0|2] =D-2{D-^) = ^-D, (5) 

with the correlation length exponent ly — 0.71 for the 
3d 0(3) transition. The 71 term is thus also irrelevant. 
However, the a term is relevant at the 3d 0(3) transi- 
tion, since it has positive scaling dimension A[a] = 0.581 
[29|. We expect this term to split the two coinciding 
0(3) transitions into two transitions, an 0(3) transition 
and an Ising transition, as was found in the Schwinger 
boson theory 16|, l30|. Again, because the 0(3) order of 
(f>i,2 implies Ising order, the latter should occur after the 
0(3) transition, i.e. at larger x. The distance in doping 
between the two transitions can be estimated by scaling, 
ignoring possible higher order singular perturbations me- 
diated by electrons: 
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FIG. 2: (a), the Brillouin zone for one-Fe unit cell at zero 
doping. The thick blue circle denotes two almost overlaping 
hole pockets |3l|, the red dashed circle represents the elec- 
tron pocket. The line which connects the electron and hole 
pockets is the (7r,0) wave vector, (b), at finite doping, the 
electron pocket expands and the hole pockets shrink, so that 
the (tt, 0) vector can no longer connect points at the Fermi 
level, (c), translation of the electron pocket by vector (tt, 0): 
If at zero doping there is a perfect overlap of the pockets, at 
infinitesimal doping there is no crossing between electron and 
hole Fermi level at all, i.e. an order with (tt, 0) wave vector 
cannot decay into particle-hole pair. 



Note that the monopoles of (j>i and <^2 are confined by 
the a term. The Berry phase for monopoles of spin-S* 
system on the square lattice is proportional to iwS, the 
monopole-composite of (j)i and 4>2 carries a trivial Berry 
phase for S* = 1 and 5 = 2 cases 
ignored hereafter. 



16, Ml, and hence is 



Note that while the 0(3) SDW order parameters 0i 
and 4>2 cannot decay into particle-hole excitations since 
the wave vector (tt, 0) does not connect pairs of points 
at the Fermi level, the same is not true for the Ising or- 
der parameter $ ~ 0i • 02 which orders at (0,0). Also 
since $ changes sign under a 7r/2 rotation and reflec- 
tion along the axis x — y, but does not break any other 
symmetry, $ couples to the two-body d— wave density 
^Q ^ EfcSign[fc2 _ fc2]cj;^^^2Cfc_^/2, and hence can de- 
cay into particle-hole excitations. The decay rate can be 
calculated using Fermi's Golden rule: 



im[x(w,q')] 



£k 



{2n\ 



2 [/(efc+g) - /(efe)](5(w - ek+q + e/c) 



X |(fc|$,|fc + g)p ~ Co-. 



(7) 



The standard Hertz-Millis [33| formalism leads to a z = 3 
theory with Lagrangian 



L^^-ni 



\UJ\ 



coq 



ciq- 



r)$„ 



(8) 



The ellipses stand for all the quartic and higher order 
terms of $ which are irrelevant at this Gaussian fixed 
point described by ([5]). Quadratic terms with singular 



factor up' jq^ or higher may occur in the expansion, but 
since the theory has z = 3 these terms are irrelevant. 
The z = 3 critical field theory was also obtained for the 
electronic nematic phase with an order parameter sim- 
ilar to $ |33|]. The critical exponents can be extracted 
directly from the field theory ([S]). For instance, in the 
quantum critical region, the specific heat and the criti- 
cal temperature of the finite temperature Ising transition 
scale as 






rpdj Z rpljZ 



\zl{d-1+z) 



Xcl 



(9) 



with d = 2. The weak interlayer coupling w^n^n+i will 
finally drive the scaling back to three dimensional be- 
havior with w ~ Jj_/J2j but its role is not considerable 
unless the 2d correlation length is long enough, i.e., if 
we are close enough to the quantum critical point. The 
spatial scaling dimension of w is A[w] = 2 at the 2d 
critical point described by Eq. [H therefore w becomes 
nonperturbative when {£,/a) '-"'' ^ 1/w, i.e. 



Xcl 



^l/(.A[H) ^ ^. 



(10) 



Within this small window, the critical scaling becomes 

(11) 






rpd/z rp 



The 0(3) order parameter 0, which can be taken as 0i, 
cannot decay into particle-hole pairs, assuming the (tt, 0) 
wave vector does not connect two points at the Fermi 
level. The Gaussian part of the Lagrangian describing 
the 0(3) transition at Xc2 has dynamical exponent z ~ 1: 



L ^ (f^q{uj^ + q^)$q + L' . 



(12) 



L' consists of quartic and higher order terms. If the quar- 
tic terms have no singularity in momentum and frequency 
space, the Lagrangian ([12]) describes a 3d 0(3) transi- 
tion. Berry phases of monopoles in this case are trivial 
for spin-1,2 16J, |30| and so are not noted. However, the 
quartic terms of the effective action may include singular 
terms like 



L' = l2\$t, 



M 



(13) 



This term can be viewed as describing the decay of |0p, 
which couples to the zero momentum charge density. 
From naive power-counting, 72 has the same scaling di- 
mension as all the other quartic terms without singular- 
ities. However, since it mixes the |(/)p field at distinct 
spatial points, the anomalous dimensions will be con- 
tributed by the two different points separately. There- 
fore its scaling dimension can be evaluated as A[72] = 
_D — 2 (_D — i) = ^ — Z?, which is again irrelevant at 
the 3d 0(3) transition. If no other more relevant quartic 



terms are present, the quantum phase transition of (p at 
Xc2 belongs to the 3d 0(3) universality class, cf. Fig. [1] 
But a thorough analysis of the quartic terms is required 
to draw a firm conclusion. 

Quantum critical points play an important role in 
transport because the electrons can scatter off the critical 
modes. We expect the Ising critical modes to contribute 
the dominant part to the low temperature resistivity, be- 
cause of its 2: = 3 soft modes and the ensuing larger den- 
sity of states at low energy. At low temperature where 
the scattering is dominated by small angle forward scat- 
tering, the resistivity is expected to scale as p ~ T^^^. 
The more general formula for the resistivity for a z = 3 
theory with Lagrangian ([5]) reads p ^ t('^+2)/^^ which 
is consistent with the well-known T^/^ law of the resis- 
tivity at the quantum critical point of three dimensional 
itinerant ferromagnetic order [34 ] . 

In summary, we have studied the SDW at (tt, 0) ob- 
served experimentally in LaOi-xFxFeAs, and its phase 
transitions. While raising the thermal and quantum fluc- 
tuations, the SDW is predicted to cede to a state with 
restored SU(2) invariance, but retaining a broken Ising 
symmetry which drives a lattice distortion. This is fol- 
lowed by an Ising transition at higher temperature or 
larger doping. The nature and universality classes of 
these transitions, and various critical exponents are dis- 
cussed. 

Note added: Fang et al. [35| have also applied thermal 
fluctuations of the J1-J2 model to the iron-based super- 
conductors. 
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